Very large stochastic resonance gains in finite sets of interacting identical subsystems 

driven by subthreshold rectangular pulses 
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We study the phenomenon of nonlinear stochastic resonance (SR) in a complex noisy system 
formed by a finite number of interacting subunits driven by rectangular pulsed time periodic forces. 
We find that very large SR gains are obtained for subthreshold driving forces with frequencies much 
larger than the values observed in simpler one-dimensional systems. These effects are explained 
using simple considerations. 
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The phenomenon of stochastic resonance (SR) seems 
to be important in a wide variety of contexts in physics, 
chemistry, and the life sciences A lot of work has 
been devoted to the study of SR both in simple Q and 
complex systems H, such as in ion channel assemblies 
or globally coupled networks of noisy neural elements [5| , 
to name a few examples. In this work, we consider a com- 
plex system formed by a finite number of N coupled noisy 
bistable subsystems, the attention to finite sets being in- 
spired by the fact that certain processes in neuroscience 
seem to involve a rather small number of subsystems Q . 

The signal-to-noise ratio (SNR) and the SR gain are 
two common quantifiers used to characterize the SR re- 
sponse of noisy systems driven by time-periodic forces. 
We will study SR effects on a collective variable of finite 
sets of interacting subunits driven by periodic rectangu- 
lar pulses. Our results show a tremendous enhancement 
of SR effects in this collective variable with respect to 
those observed in single unit systems. 

Let us consider a set of N interacting subsystems, each 
one of them characterized by a single degree of freedom 
Xi (i — 1, . . . , N), whose dynamics is governed by the 
Langevin equations @, H, Q 
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where &(t) are Gaussian white noises with zero average 
and (£j(i)£j(s)) = 2DSij5(t — s), 9 is the parameter defin- 
ing the strength of the interaction between subsystems, 
and F(t) is an external driving force of period T. In this 
work, we will restrict ourselves to forces of the type [ijj 



F(t) = 



A ; < t < t c , 

; t c < t < T/2, 

-A ; T/2<t< T/2 + t c , 

; T/2 + t c < t < T. 



(2) 



The parameter r — 2t c /T, usually called duty cycle, mea- 
sures the fraction of a period during which this driving 
force has a nonvanishing value. In addition, we will only 
consider subthreshold amplitudes A, so that the driving 
force ((2]) cannot induce sustained oscillations between dy- 
namical attractors in the absence of noise. The model 



described by Eq. JT]) (without the external periodic driv- 
ing) was used years ago by Kometani and Shimizu 0] as 
an empirical model to describe muscle contraction. Later 
on, Desai and Zwanzig [§[ gave a more detailed statistical 
mechanical description in the asymptotic N — > oo limit 
and used it to model order-disorder transitions. The ad- 
dition of an external driving can, in principle, be used 
to describe the phenomenology of a forced contracting 
muscle. 

We focus on the collective variable S(t) defined as 
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S(*) = ^5>(*), 



(3) 



which has previously been used Q in the global analysis 
of coupled bistable systems. It might be considered as 
the total output process of a parallel array of N inden- 
tical interacting subunits, subject to independent noise 
sources, £j(t) and the same external forcing F(t). Its 
signal-to-noise ratio (i? ut) is defined in the usual way 
as, 
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dr C co h(T) cos(fir) 
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Ql = - drCi nco h(r)cos(Qr). 
i" Jo 



(4) 
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where Q = 2vr/T, C coh (r) = i / Q r dt(S{t + r)} ^ (S '(t)> oo 
and C inco h(r) = C(r) - C co h(T) with C(r) = 
i r T ^(5(t + r)^)) oo . 

For a set of N coupled linear oscillators driven by the 
external driving force F(t) and subject to the noise terms 
£i(t) as in Eq. (pj) , the SNR of the corresponding col- 
lective process, R^h coincides with that of the random 
process formed by the arithmetic mean of the individ- 
ual noise terms plus the deterministic driving force 
F(t), namely, F(t) + £(f) with £(t) - N' 1 £]Li &(t). 




FIG. 1: Dependence of the signal-to-noise ratio R ou t and the 
SR gain G for a set of N = 10 identical subsystems and an 
external driving of frequency £1 = 0.01 for several values of the 
duty cycle r: 0.1 (open circles), 0.2 (crosses), 0.3 (triangles), 
0.4 (squares), and 1 (stars). Inset shows the SR gain as a 
function of r for a fixed noise strength D — 0.08. The lines 
have been drawn as a guide to the eye. 



The process £(t) is a Gaussian white noise of effective 
strength D/N. Then, it is easy to prove that 



{L) _ 2A 2 N[l - cos(?rr)] 
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Thus, for our nonlinear case, it seems convenient to ana- 
lyze the SR gain, G, defined as 0] 



G = 



Re 



R 



(L) 



(8) 



which compares the SNR of a non-linear system with 
that of a linear system subject to the same stochastic 
and deterministic forces. 

We have carried out extensive simulations of the 
Langevin equations §Q with the external driving |2|). 
In all cases reported here the coupling strength is fixed 
to 9 = 0.5 and the subthreshold driving amplitude to 
A = 0.3. There is nothing special about this particular 9 
value. The qualitative results would be the same for any 
other value of 9 ^ 0. 

In Fig. Q] we show the collective signal-to-noise ratio 
and SR gain as a function of the noise strength D for a set 
of N = 10 identical subsystems, a driving fundamental 
frequency £1 = 0.01 and several values of the duty cycle 
r. It can be seen that, while the SNR curves are nearly 
identical for r < 0.4, the SR gain increases drastically, 
reaching a very large value for r — 0.1. This dramatic 
increase is easily understood by taking into account the 
observed almost constant behavior of i? ou t and the fact 
that i?out j Eq. ([7]), decreases monotonically with r. How- 
ever, for sufficiently small r, i? ut must decrease faster 
than because the interval t c becomes smaller than 

the time it takes for the system to react to a constant 
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FIG. 2: Dependence of the denominator and the numerator 
of the collective signal-to-noise ratio for the same parameter 
values as in Fig [1] The lines have been drawn as a guide to 
the eye. 



force of amplitude A, and thus, the driving produces al- 
most no effect in the system. This behavior is shown in 
the inset of Fig. [T] where it can be seen that the SR gain 
decreases for r < 0.08. Also, as seen in Fig. [TJ when 
r = 1 (rectangular driving signal) , G reaches a peak at a 

noise value D ps 0.2. Even though R^l i s as large as it 
can be, the huge increase of i? ut is enough to overcome 
the increase in i?^, yielding a substantial value for the 
SR gain. 

The very large gain values in Fig. [T] for pulses with 
short duty cycles are observed only for a small range of 
noise strengths around D « 0.08. In Fig. [2] we present 
the behavior of Qi and Q u with D. Notice that around 
D w 0.08 there is a strong reduction of two orders of mag- 
nitude in the level of fluctuations of the collective variable 
as measured by Qi. Therefore, the large SR effects quan- 
tified by the SNR and the SR gains are essentially due 
to the very large reduction of the fluctuation spectrum of 
the output signal at the fundamental driving frequency 
for a range of noise values. 

We have also analyzed the SNR and the SR gain for 
different values of f2. Our results for i? out and G as a 
function of D are depicted in Fig. [3] for a set of N = 
10 interacting identical subunits driven by rectangular- 
pulses with duty cycle r = 0.1 and several values of the 
fundamental driving frequency Q: 0.015 (open circles), 
0.02 (crosses), 0.03 (triangles), and 0.04 (squares). As we 
increase the driving frequency O, the SR gain is gradually 
reduced. For a sufficiently large driving frequency (pulses 
of very short duration) , gains larger than unity are not 
observed for any value of the noise strength. 

As mentioned above, the large SNR values observed in 
Fig. [2 are related to the existence of a sharp minimum of 
Qi at a certain noise value. As observed in Fig. [2l this 
value is D w 0.08 for a driving signal with short duty 
cycle, amplitude A = 0.3 and fundamental frequency 

= 0.01. For other amplitude and frequency, the lo- 
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FIG. 3: Dependence of the signal-to-noise ratio and the SR 
gain for a set of N — 10 identical subsystems and duty cy- 
cle r = 0.1 for several values of the driving frequency Q: 
0.015 (open circles), 0.02 (crosses), 0.03 (triangles), and 0.04 
(squares). Solid lines have been drawn as a guide to the eye. 
The horizontal dashed line marks the unity for the SR gain. 



cation of the minimum will change, although the mech- 
anisms leading to the existence of this minimum will be 
qualitatively the same. 

To understand the Qi dependence onD, it is important 
to analyze in detail the dynamics imposed by the external 
driving force in Eq. ([2]). For the D and 9 values of inter- 
est to the present discussion, simulations show that when 
the driver is absent, the collective variable S{t) performs 
a noise induced random movement between the two sym- 
metric attractors of the dynamics, which we will regard 
as located at ±Sq. During a time interval of duration t c , 
the external force F(t) = A favors the positive attrac- 
tor. Thus, if S(t) was at the negative attractor at the 
beginning of the interval, the external forcing will drive 
it to the positive one in a random time that we will de- 
note by Ti. Certainly, different realizations of the noise 
will yield different values of Ti. Running simulations 
with many independent trajectories, we have computed 
the probability Prob(Ti < r) = /(r) that the variable 
S(t) has jumped before a time r to the attractor favored 
by a constant driving of amplitude A. In Fig. 2] (a) we 
present /(r) for the set of parameters relevant to the dis- 
cussion: A = 0.3, D = 0.08, and N = 10 (solid line). 
It can be seen that for the case r = 0.1 and = 0.01, 
(thus, t c w 31.4), a transition between the attractors for 
t = t c is performed with probability almost unity. Then, 
during the rest of the half-period, an interval of duration 
T/2 — t c , the external force is zero and the system is free 
to jump between the attractors due solely to noise. Let 
us now denote by T2 the random time it takes to jump 
from one attractor to the opposite one when F(t) = 0. 
Figure |4] (b) shows the probability Prob(T 2 < r) = g(r) 
that this jump has taken place before a time r vs. r. 
Since for r = 0.1 and Q = 0.01 we have T/2 - t c « 282.7, 
it can be checked in Fig. Q] (b) that almost no tran- 
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FIG. 4: (a) Probability that the collective variable S(t) has 
jumped before a time r to the attractor favored by an exter- 
nal constant driving of amplitude A = 0.3 starting from the 
opposite attractor. The solid line corresponds to a system 
with D = 0.08 and N = 10, the dotted line to D = 0.08 and 
N = 1, and the dashed line to D = 0.04 and N = 10. (b) 
Probability that S(t) jumps before a time r to the opposite 
attractor in the absence of driving. The solid line corresponds 
to a system with D = 0.08 and TV = 10 (being almost zero for 
the range of times shown), the dotted line to D = 0.08 and 
N = 1, and the dashed line to D = 0.16 and N = 10. 



sitions take place during this time interval under these 
conditions. We can carry out the same analysis for the 
symmetric situation during the second half-period of the 
driving force. Consequently, S(t) performs a neat tra- 
jectory between its attractors with transitions induced 
systematically every half-period by the external driving 
when it has a nonvanishing value. As a result, we would 
expect Ci nco h(0), itself an average of the second cumulant 
of S(t) over a period, to be of the order of the effective 
noise D/N, and Ci n coh('?") short-lived. This is, in fact, 
what is observed in Fig. [5] (see inset) . 

The curves in Fig. HJa) indicate that the probabil- 
ity of transitions between the attractors before t c out of 
the most unstable attractor is smaller for D = 0.04 than 
for D = 0.08. This is to be expected as the probabil- 
ity of transitions out of an attractor goes to zero when 
D — > with subthreshold driving amplitudes. As a re- 
sult, Cincoh(O) and the decay time of Ci nco h(T) increase 
as D gets smaller than D = 0.08. Consequently, notic- 
ing its definition in Eq. Qi has to increase as D is 
reduced from D = 0.08. 

Fig.[5]shows that for noise values larger than D w 0.08 
the Qi dependence on D for inputs with short duty cy- 
cles is rather different from the one observed for a rect- 
angular input (r = 1). This fact indicates that for short 
duty cycles and these larger noise values, the time inter- 
vals during which F(t) — are crucial for the system 
response. By comparing the g(r) curves for D = 0.08 
and D = 0.16 in Fig. Q] (b), we see that the probability 
of transitions between attractors during the time interval 
T/2 — t c is negligible, whereas there is a non- negligible 
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FIG. 5: Temporal behavior of the incoherent part of the 
correlation function of the collective variable for a system of 
N = 10 identical subsystems with D — 0.08 and r = 0.1. 
The solid lines correspond to £1 = 0.01 and the dashed line to 
f2 = 0.03. The inset is a magnification of the case Q — 0.01 
displayed for the sake of clarity. 



probability for D = 0.16. Therefore, the bimodal charac- 
ter of the probability distribution for S(t) is enhanced as 
D increases from D w 0.08 to D w 0.16 and one can then 
conclude that Ci nC oh(0) and the decay time of Ci nco h(r) 
increase as D increases. Consequently, for short duty cy- 
cles, Qi also increases as D increases within the interval 
D « 0.08 to D « 0.16. On the other hand, for r = 1 and 
this same range of noise values, the situation is different 
because the driving signal never vanishes. Then, tran- 
sitions out of the most unstable attractor happen more 
frequently as noise increases [/(r) saturates at 1 earlier], 
and this leads to a decrease of both Cmcoh(O) and the 
decay time of Ci nC oh(T) as D increases, with the corre- 
sponding decrease in Qi. 

Similarly, raising the driving fundamental frequency 
has a drastic effect on the behavior of Cincoh^) with re- 
spect to its optimal behavior discussed above for D = 
0.08 and r — 0.1, as depicted in Fig. [5] for two fre- 
quency values: £1 = 0.01 and f2 = 0.03 . The picture 
that emerges for f2 = 0.03, differs considerably from that 
observed for = 0.01. For the higher frequency case, 



with t c w 10.4, the value of the solid line at r = 10.4 
in Fig. Q] (a) indicates that, for a considerable number 
of trajectories, the driving force is not able to induce a 
transition during the time interval t c . Consequently, the 
second cumulant of S(t) becomes of the order of the dis- 
tance between the attractors. On the other hand, for the 
lower frequency, transitions occur for the corresponding 
r = t c w 31.4 with almost total certainty. This leads to 
a Cincoh(O) for ft = 0.03 much larger than for = 0.01. 
Additionally, since still almost no transitions occur at the 
intervals when the driving is absent, as shown in Fig. 0] 
(b), nor obviously when the driving favors the initial at- 
tractor, any unsuccessful transition is carried on until the 
next period, and we would expect a long-lived C?i nco h(r). 
Namely, noise induced correlations persist during a few 
driving periods. These observations justify the full be- 
havior observed in Fig. [5l 

In conclusion, we have analyzed the enhancement 
of SR effects in a collective variable characterizing the 
response of finite sets of interacting noisy subsystems. 
The cooperative effect of noise and nonlinearity in finite 
sets of interacting, driven, bistable subunits reflects in a 
substantial decrease in the noise level of the collective 
output, S(t), with respect to that observed in the output 
of a single independent unit. Very large values of the 
SR gain can be achieved for trains of short rectangular 
pulses of the type defined by Eq. ([2]). Even though the 
results presented here have been obtained with external 
pulses with very brisk changes of their amplitudes at cer- 
tain instants of time, similar results can also be observed 
when the pulse amplitude changes continuously, as long 
as there are time intervals within a period where the 
amplitude changes very slowly, followed by short time 
intervals where the amplitude changes very drastically. 
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